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The kinematic identification of a thick stellar disc in M31 
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ABSTRACT 

We present the first characterization of a thick disc component in the Andromeda 
galaxy (M31) using kinematic data from the DEIMOS multi-object spectrograph in- 
strument on Keck II. Using 21 fields in the South West of the galaxy, we measure 
the lag of this component with respect to the thin disc, as well as the dispersion, 
metallicity and scale length of the component. We find an average lag between the 
two components of (Av) — 46.0 ± 3.9 kins" 1 . The velocity dispersion of the thick 
disc is (Tthick = 50.8 ± 1.9 kms -1 , greater than the value of dispersion we determine 
for the thin disc, athin = 35.7 ± 1.0 km. B . The thick disc is more metal poor than 
the thin disc, with [Fe/H] spec = —1.0 ± 0.1 compared to [Fe/H] spec = —0.7 ± 0.05 
for the thin disc. We measure a radial scale length of the thin and thick discs of 
h r = 7.3 ± 1.0 kpc and h r — 8.0 ± 1.2 kpc. From this, we infer scale heights for both 
discs of 1.1±0.2 kpc and 2.8±0.6 kpc, both of which are ~2-3 times larger than those 
observed in the Milky Way. We estimate a mass range for the thick disc component 
of 2.4xlO lo M < M*_thick < 4.1 x 10 10 M Q . This value provides a useful constraint 
on possible formation mechanisms, as any proposed method for forming a thick disc 
must be able to heat (or deposit) at least this amount of material. 



1 INTRODUCTION 



Roughly 70% of bright ga laxies observed at re d shift z = 
possess stellar d i scs (e .g. iHammer et al. 1 120051; iPark et al.l 
120071 : IChoi et alj|2007l and iDelgado-Serrano et all |201Cl 11. 
including our own Galaxy and its two largest neighbours, 
M31 and M33. From this, we can infer that spiral mor- 
phologies are the dominant configuration for galaxies 
viewed at the present epoch. Under the formalism of 
hierarchical structure formation, galaxies are believed to 
evolve into their present forms via the accretion of, and 
mergers with, smaller systems. The effect of this process 



1 The data presented herein were obtained at the W.M. Keck 
Observatory, which is operated as a scientific partnership among 
the California Institute of Technology, the University of Califor- 
nia and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support 
of the W.M. Keck Foundation. 

2 Based on observations obtained with MegaCam, a joint project 
of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Tele- 
scope (CFHT) which is operated by the National Research Coun- 
cil (NRC) of Canada, the Institute National des Sciences de 
l'Univers of the Centre National de la Recherche Scientifique of 
France, and the University of Hawaii. 



on the seemingly dynamically stable stellar discs we see in 
Milky Way (MW) type galaxies is still largely uncertain. 
The ability of these fragile objects to survive a "major 
merger" event (i.e. a merger with a system > 1/3 of the 
host's mass) is something that is still debated. Such major 
mergers are thought to be cosmologically common, with 
~ 70% of all galaxies with a halo mass of M ~ 10 12 M 
having experienc ed at least one major interaction withi n 
the past 8 Gyrs (|Stewart et al J 12008 : [Purcell et al.1 120091 '). 
Thus it has been argued that galaxies that possess thin 
stellar discs at z = could not have experienced a major 
merger within the last 10 Gyr wi t hout the disc being 
destabilized jToth fc Ostrikerl Il992l: IWalker et al.1 1 19961 : 
IStewart et"afl l200Sl : iPurcell et alj 120091 1 . This poses a 
significant challenge to our understanding of the formation 
of disc galaxies like the MW and M31. Recently, several 
authors have argued that these thin discs could survive 
such an event if the me r ging system i s sufficiently gas rich 
jRobertson et al.l [200(1: iBrook et al.l 120071; iHopkins et al.1 

although 



2009; 



Stewart et al.l 



2009: 



Brooks et al 



2009), 



the disc would still undergo heating, resulting in a thicker 
disc than that observed presently in the MW. In addition 
to the effect of major mergers on the structure of discs, 
galaxies viewed at the present epoch have undergone (and 
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are still undergoing) many smaller "minor" mergers which 
are not sufficiently massive to destroy thin stellar discs, but 
are thought to kinematically heat them, causing them to 
flare outwards and create a second, t h ick disc component 
llQuinn et all 1 19931; iRobin et al l Il99rj; IWaiker et al.l Il996l; 
Velazquez fc White! Il999l; IChen et all l200ll; ISales et all 



Elmegreen fc Elmegreen! 120061 ; lYoachim fc Dalcantonll20(56l . 



2009; Vill alobos fc H clmi 20091; iPurcell et al.l l2010h . Other 



physical processes are also thought to heat up and 
thicken the thin disc, including the ac cretion of a satel- 
lite on a radial orbit about its host (|Abadi et all 120031 ; 
I Read et al.l 12008ft . internal heating within the disc from 



dVillumsen 19851; Carlbera 


19871; Sellwood & Binncv 


2002; Hanninen & Flvnnl 


20021; Benson et alj 


2004; 


Havashi & Chiba 20061: Kazantzidis et alj 


2006; 


Roskar et al. 20081; ISchonrich & Binnevl 12009a; 


Loebman et al.l 12010ft. Thick 


discs may also have formed 



thick, with significant star formation occurring above 
the mid-pla ne of the galaxy o r with large ini tial velocity 
dispersio ns dBrook et al] |2004| ; iKroupal 120021 ) . In recent 



work by Ros kar et alj (|2010h . they suggest that in-situ 
formation could also occur if the stellar disc is misaligned 
with the hot, gaseous halo. This misalignment results in 
a significant warping of the outer disc, and subsequent 
star formation within this warp results in a low metallicity 
thick disc. Finally, it is also possible that a number of these 
mechanism will act in conjunction. In particular, it has been 
suggested by a number of authors that secular growth from 
internal heating may be significantly enha nced by minor 
merger events via swing a mplification (e.g. ISellwood et al.l 
Il998l : iDubinski et aTll200Sft . as these processes often occur 
simultaneously. As such, it makes little sense to treat these 
two scenarios as separate processes. 

With so many potential mechanisms capable of pro- 
ducing thickened stellar discs, just how common are 
thick discs in spiral galax ies at the present epoch? 
iDalcanton fc Bernstein! (|2002l ) claim that thick disc forma- 
tion is a universal feature of disc formation, and as such 
should be observed in all spiral galaxies. Whether such discs 
are formed predominantly via one mechanism, or a mixture 
of them is still uncertain, and disentangling the various for- 
mation scenarios from one another in present data sets has 
proven difficult. 

In the MW, the existence of a thic k disc has long 
been known, and was first identified by iGilmore fc Reidl 
(1983). Subsequent spectroscopic studies of this compo- 
nent have shown it to be kinematically distinct from 
the thin stellar disc, with the thic k disc lagging be - 
hind the thin disc by ~ 50kms _1 (|Carollo et alj 120101 ) 
and having a larger velocity dispersion than the thin 
disc. This thick component also seems t o be composed 
of older, more metal deficient stars (e.g. IChiba fc Beers! 
2000; IWvse et al.l 12006ft . However the observed properties 
of the thick disc, such as scale height, length and velocity 
dispersion, tend to vary dep ending on the survey sample 
and t racer population used (lJuric et al. 20081 ; llvezic et al.l 
|200S| ; ICarollo et ai]|2010l ; Ide Jong et alj|2010ft . As such, the 
origin of the MW thick disc is still a subject of great 
debate in the literature. Thick discs have also been ob- 
served in a number of edge on spiral galaxies (e.g. iBursteinl 
19791 ; lTsikoudlll979l; Ivan der Kruhj|l984 IShaw fc Gilmorel 
~ ~E>okkum et alj 1994 IDalcanton fc Bernsteinll2002l ; 



l2008bH al). and spectroscopic observations of these objects 
also show the thick discs to be composed of older stars than 
their corresponding thin discs. However, as these galaxies 
are all located at distances greater than ~10 Mpc from the 
MW, one cannot obtain spectra for individual stars, and 
must instead rely on the integrated spectral properties of 
RGB stars. Obtaining spectra with a high enough signal-to- 
noise (S:N) to discern velocity dispersion profiles and reli- 
able metallicities is also challenging, making it impossible 
to distinguish between the various formation mechanisms 
for these structures. 

If thick stellar discs are universal amongst spiral galax- 
ies, and are formed by mergers with, or accretions of, 
satellites, one might expect to see such a structure in 
M31. This neighbouring galaxy is considered to be a "typ- 
ical" spiral galaxy when compared with other local exter- 
nal disc galaxies (|Hammer et al 1 l2007h . It is thought to 
have had an active merger history, and a recent panoramic 
photometric survey by the Pan-Andr omeda Archaeolog- 
ical Survey collaboration (PAndAS, iMcConnachie et al.l 
2009) has shown the halo of this galaxy to be littered 
with tidal streams from interactions with in-falling satel- 
lites. These incl u de the Giant Sou thern Stream (GSS, 
llbata et all l200ll; iGilbert et all 120091') and streams A, B 



1989; 



C, D and E llbata et al.1 120071 ; IChapman et al.l 120081 ; 
IMcConnac hie et al. 2009 ). The outer disc of M31 is ver y per- 
turbed dFerguson et al.ll2002l ; [Richardson et al.l 12008ft . sug- 
gestive of some tidal interaction. M31 is also host to 25 
known dSph and 4 dE satellites, at least 2 of which (NGC 
205 and M32) show evidence for significant tidal interac- 
tion (IChoi et al.ll2002l; IMcConnachie et al.|[2004l ; lGeha et all 
120061 ; iHowlev et al.1 12008! '). Therefore the possibility of nu- 
merous interactions between the disc of M31 and its satellite 
population seems highly likely. IMcConnachie et alj (|2009ft 
also present evidence for a recent interaction between 
M31 and its neighbouring spiral galaxy, M33, which could 
have significantly distorted and heated the M31 disc, giv- 
ing rise to a thick disc component or substantial sub- 
structure in the outer disc. Other groups have postulated 
links between t he formation of bulg e s and thick discs in 
spiral galaxies dMelendez et al.l 12008); [Hopkins et al. 1 120081 ; 
iBournaud et al.ll2009l ; iBensbv et allboicft . and as M31 is 
know n to have a reasonably massive bulge $sEg± 

ia et al.l 

l2010ft . it is an interesting candidate for hosting a thick disc. 
Despi te its high inclination to us along the line of sight 
(77°, IWalterbos fc KennicuttJ ll988T ). M31 is not seen suf- 
ficiently edge-on to allow us to look for such a population 
using photometry. Therefore to look for evidence of a thick 
disc in M31, we must search for it via its kinematic signa- 
ture, using spectroscopy. Given its proximity to us (785 kpc, 
IMcConnachie et al.ll2005ft . M31 is an ideal target for spec- 
troscopic observations as we are able to resolve and obtain 
reliable velocities for individual Red Giant Branch (RGB) 
stars, and it has an advantage over our own galaxy as we 
are afforded a panoramic view, whereas in the MW we are 
hampered along various lines of sights by confusion from the 
disc and the bulge. 

Since 2002 our group has been conducting a system- 
atic spectroscopic survey of M31, including the disc, halo 
and regions of substructure using the DEIMOS instru- 
ment mounted on the Naysmyth focus of Keck II (105, 
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Figure 1. Map showing the location of fields within our M31 
survey. Fields selected for study in this work are labelled and 
highlighted in red. The outer ellipse shows a segment of a 55 kpc 
radius ellipse flattened to c/a= 0.6, and the major and minor 
axis are indicated with straight lines out to this ellipse. The inner 
ellipse corresponds to a disk of radius 2°, (27 kpc), with the same 
inclination as the main M31 disk. 



Chapman et~aT]|2005l , l200fj , 120071 . l200St ICollins et al1l2009l . 
2010T ). The data from this survey gives us an ideal opportu- 



nity to identify a thick disc if present. In fact, in their study 
of M31's extended disc using this same data set, 105 iden- 
tified a population lagging behind the thin disc which they 
excluded from t heir study that they te rmed a 'thick disc-like' 
population and I Chapman et al.l (|2006t ) briefly examined this 
component as a function of radius but were unable to com- 
ment on its global properties. In this work, we discuss the 
results from an in depth study of this population, analysing 
its kinematics and chemistry, comparing them to M31's thin 
and extended discs, the thin and thick discs in the MW, and 
those seen in other galaxies, and comment on possible for- 
mation scenarios for this component. The paper is set out 
as follows; in §2 we discuss the known structure of M31, §3 
focuses on our spectroscopic survey of M31 and discusses 
field selection and analysis techniques. We present our re- 
sults in §4 and discuss their implications in §5. We conclude 
our findings in §6. 



2 THE BULGE, DISCS AND HALO OF M31 

The first recorded observation of M31 was made in 964 AD 
in the 'Book of constellations and fixed stars' written by the 
Persian astronomer, Abd al-Rahman al-Sufi, who described 
it as 'a small cloud' in the night sky. In the centuries that 
have passed since, M31 has been a popular target for as- 
tronomers, and much has been learned about its structure. 
M31 is a spiral galaxy of SA(s)b type, with a significant 



bulge, a classical thin stellar disc, a vast extended stellar 
disc and a metal poor halo. In this section, we outline the 
properties of each of these components. 

First we discuss the bulge component. Numerous stud- 
ies have show that M31 possesses a classical bulge, with 
a Sersic index of ~ 2 and an effective radiu s of 1.93 kpc 
|Kormendv fc Benderlll999l ; ISeigar et al.ll2008T ). It is largely 
supported by rand om motions, although recent work by 
ISaglia et all |2010h has found evidence for rotation in the 
innermost regions. Sagiia et al. also find the bulge to be 
dominated by an old stellar population (age > 12 Gyrs) of 
roughly solar metallicity, with a large velocity dispersion of 
166-170 km s -1 . This component is dominant out to about 
8' (~ 2 kpc) , at which point the disc begi ns to dominates th e 
surface brightness pr ofile of the gal a xy ((Saglia et aDl2010T ). 
however according to lMerrett et al. the bulge can be 

traced out as far as 10 effective radii, equivalent to ~ 15 kpc 
meaning that some of our innermost disc fields may be sub- 
jected to minor contamination from this component. 

Studies of the thin stellar disc of M31 have been per- 
form e d by a number of auth o rs (e.g . IWalterbos fc Kennicutia 
19881: iFereuson fc Johnson! l200ll" 105: llbata et al.1 120071 



McConnachie et all l2009h . and have challenged our pre- 
vious notions of the struc ture of classical discs. With 
a scale length of 5.9 kpc (jWalterbos fc KennicuttJ 1 19881 . 
corrected for assum ed distance to M31 of 785 kpc, and 
iMerrett et al] [2006), it is more extensive than that of our 
Galaxy, and also appears t o be forming stars at a lower rate 
[jWalterbos fc BraurJfl994T ). And while it is a characteristic 
feature of the surface brightness pr ofiles of stellar discs to 
steep l y decline at 3-4 sca le lengths l)van der Kruit fc Searld 
Il98ll : [Pohlen et ai]l200d ). which corresponds to 18 -24 kpc 
in M31, a spectroscopic study by llbata et ap ( 2005) uncov- 
ered a vast, extended disc component that can be traced 
out to distances of ~ 40 kpc (~ 8 scale lengths) from the 
centre of the galaxy that has an exponential surface den- 
sity profile that is very similar to the inner disc. While 
this structure is rather clumpy, on average it appears to 
follow on smoothly from the classical inner disc, although 
perhaps with a slightly larger scale length of 6.6±0.6 kpc 
and with a slight lag behind circular velocities at large radii 
( (Av) = 20kms _1 , 105). It is dynamically cold, with a ve- 
locity dispersion ranging from 20-40 km s _1 , leading 105 to 
conclude that it is likely not a thickened disc. Whether this 
extended component is truly separate from the classical thin 
disc, or merely an extension of it that shows some evidence 
of heating and warping at larger radii where the disc is more 
sensitive to perturbations from mergers and interactions, is 
unclear. Owing to the similarity of the thin and extended 
discs, we will refer to them both as the 'thin disc' through- 
out this paper. Where we wish to make a distinction between 
the two, we shall use the terminologies 'classical' and 'ex- 
tended' disc. 

The presence of a smooth, pressure supported metal- 
poor halo in M31 el uded detection until very recently. 
In 20 06, two groups |Chapman et al.l 120061 : iKalirai et al.l 
2006) independently identified such a component using the 
DEIMOS instrument on the Keck II telescope. Centred on 
the systemic velocity of M31, with a central velocity dis- 
persion of 152 km s" 1 , and showing no strong evidence of 
rotation, both groups found this component to be metal 
poor with an average metallicity of [Fe/H] = —1.4 ± 0.2 
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l|Chapman et all 120061 ). iKalirai et all (2006) were able to 
trace this component out as far as 165 kpc from the cen- 
tre of the galaxy although there is an inevitab le confusion 
with the halo of M33 at these large distances (Ibata et al.l 
120071 ; iKoch et al.ll2008l ; IMcConnachie et al.ll200gF r 

The halo of M31 is also a known host to a number of 
kinematic substructures, such as the GSS, the tangential 
streams that cross the SE minor axis, the western shelf and 
a wealth of substructure in the NE of the galaxy that is 
thought to be linked to the GSS. In the following analysis, we 
will carefully consider the kinematics of all these components 
to ensure any thick disc sample that we define is free from 
contamination by any of these sources. We shall discuss this 
in greater detail in §3. 



3 OBSERVATIONS AND FIELD SELECTION 

A detailed description of the observational methodology and 
target selection employed in the survey is given in 105, 
which we briefly summarize here. Using Colour-Magnitude 
Diagrams (CMDs) from both the Canada France Hawaii 
(CFHT) and Isaac Newton (INT) telescopes ([Ferguson et all 
120021 ; llbata et all 120071 ; IMcConnachie et al]|2009l ), we se- 
lected targets for observation by prioritising Red Giant 
Branch (RGB) stars in M31 with 20.5 < i < 21.2 and colours 
1.0 < (V — i)o < 4.0 (priority A), then filling the remainder 
of the masks with stars with I < 22.0 that are unsaturated 
(priority B), where the V and I colours are transformed from 
their native g and i colour s usi ng the relations de scribed in 
IMcConnachie et all l|2004l ) and llbata et all (|2007h . We used 
a combination of standard DEIMOS multislit mode for low 
density fields, such as the halo, and our own minislitlet ap- 
proach which allowed us to target > 600 stars per mask in 
more crowded regions, such as the disc. Our observational 
setup covers the range of the Calcium Triplet (Ca II) lines 

o 

at 8498, 8542 and 8662 A, a prominent absorption feature 
that can be used both to measure radial velocities, and as a 
metallicity indicator. To obtain velocities, we cross-correlate 
all observed stars with a template Ca II spectrum. We esti- 
m ate the errors on our v elocit ies by following the procedures 



m ate t he errors on our v elocities by following the procedures 
of ISimon fc Gehal (|2007l ) and lKalirai et all l|2010l ). First, we 
make an estimate of our velocity uncertainties for each ob- 
served star using a Monte Carlo method, whereby noise is 
randomly added to each pixel in the spectrum, assuming a 
Poisson distribution for the noise and the velocity is recal- 
culated using the same cross correlation technique described 
above. This procedure is repeated 1000 times, and then the 
error is calculated to be the square root of the variance of 
the resulting mean velocity. We combine this error with a 
systematic error, e, which contains information on any errors 
we may not have accounted for (for example, wavelength cal- 
ibration error, misalignment of the mask etc.). For the fields 
observed with the 600 line/mm grating, we evaluate this er- 
ror directly by using repeat measurements in fields 231Dis 
and 232Dis, a total of 332 stars. We define the normalised 



error, crjv as: 



(JJV 



V 1 — V2 



(a 1 2 +a| + 26 2 ) 1 /2 



(1) 




ct n (km/s) 

Figure 2. A histogram showing the normalized error distribution 
for repeat measurements of the same stars in two of our fields that 
were observed with the 600 lines/mm grating. The normalized 
error, (Xjv incorporates the velocity differences between the repeat 
measurements (vi and V2), and the Monte-Carlo uncertainties 
calculated for each observation (cti and 0-2). In order to reproduce 
a unit Gaussian distribution for our uncertainties, we also include 
an additional error term, e, which accounts for any systematic 
uncertainties we have not included. We find e=5.6 kms -1 for this 
setup. For fields usin g the 1200 lines/mm setup, we using the 
ISimon fc Gehal (§007) value of e = 2.2 km s" 1 . 



required in order to reproduce a unit Gaussian distribution 
with our data (shown in Fig. [2} This gives us a systematic 
error for this setup of e = 5.6kms~ , slightly lower than the 
value of e = 6.2 km s -1 derived bv lCollins et al.l (|2010l ) for 
the same setup, though we note that their measurement was 
based on repeat observations of 47 stars, compared with our 
much larger data set of 332 stars. The typical uncertainties 
for these measurements above a threshold of S:N = 3, are 
5-10 kms" 1 . 

The Ca II features also provide us with a method for 
measuring the spectroscopic met allicity of our ob s erved sam- 
ple. Following the pro cedure of iRutledge et al] l|l997h and 
(Battarii a~etaLl |200Sl ). we fit Gaussian functions to the 
three Ca II peaks to estimate their equivalent widths (EWs), 
and calculate [Fe/H] using equation (1) 



[Fe/H] = -2.66 + 0.42[ECa + 0.64(V RG b - V HB )] 



(2) 



where Vi and 1)2, o\ and 02 are the velocities and errors of 
each measurement pair, and e is the additional random error 



where ECa=0.5EW 8 498+EW 8 542 +0.6EW8662, Vrgb is the 
magnitude (or, if using a composite spectrum, the aver- 
age magnitude) of the RGB star, and Vhb is the mean V- 
magnitude of the horizontal branch (HB). Using Vhb — Vrgb 
removes any strong dependence on distance or reddening 
in our calculated value of [Fe/H], and gives us the Ca II 
line strength at th e level of the HB. F or M31, we set this 
value to be 25.17 l|Holland et al.lll996l ). We note that this 
assumed value is sensitive to age and metallicity effects, see 
[Ch en et al. I l2009l for a discussion, however owing to the large 
distance of M31, small differences in this value within the 
disc of M31 will have a negligable effect on metallicity calcu- 
lations. For individual stars, these measurements carry large 
errors (> 0.4 dex), but the errors are significantly reduced 
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Figure 3. Here we show a number of HI rotation c urves for M31. 
Throughout this work, we use results based on the lChemin et all 
(120091 1 work, shown as fille d black circles. We also show rotation 
curves from ICorbelli et al. (2010) and 105, which differ slightly 
from the Chemin et al. 2009 curve i n the outermo s t regi ons. We 
find that using these curves vs. the IChemin et alj l l2009h results 
do not affect our results. 



when stacking the spectra into a composite in order to mea- 
sure an average metallicity for a given population. 



3.1 Field selection and sample definition 

In order to detect a thick disc component in M31 kine- 
matically, we need to measure the velocities of stars within 
our sample relative to some model for the velocity of stars 
within the thin disc of the galaxy. If a thick disc is present, 
we should observe a population that lags behind the thin 
disc in terms of its rotational velocity. The component is 
also expected to have a larger velocity dispersion than the 
thin disc. This is observed in the MW, where the thick disc 
lags the thin by betw een 20-50 kms" 1 l|Chiba fc Beersll2000l: 
ISoubiran et~afll2003h and has an average rotationa l velocity 
dispersion of ov^=57 kms -1 (|Carollo et alJl2~01Ch . For the 
purposes of this work, we shall use an updated version to 
the disc model of 105. In this model, we assume circular or- 
bits for all stars about the centre of M31 a nd interpolate 
their velocities from the HI rotation curve of[Ch emin et al.l 
(2009), which is shown in Fig. [3] as the solid black points. 
This rotation curve differs from that adopted by 105, par- 
ticularly in the oute rmost regions . The y used a compila- 
tion of CO data f rom Klypin et al.l (|2002t ) and HI data from 
iBrinks fe Burtonl (|l984l ). which we also show in Fig. [3] as 
re d triangles. We also show the HI rotation c urve derived 
bvlCorbelli et all (|2010r i from a WRST survey l|Braun et all 
2009) as blue squares. Using either of these rotation curves 
as opposed to that of Chemin et al. leads to differences in 
our interpolated velocities of order a few-20 kms -1 , however 
there is a negligible effect on the dispersions within particu- 
lar populations, and so the adoption of any of these curves 
would give us consistent results when analysing the global 
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Figure 4. A contour map of the expected velocities of stars in 
circular orbits in the disc of M31. This was constructed using our 
simple model as discussed in §3. 



properties of the stellar dis cs in this work. We assume a n 
inclination for M31 of 77° l|Walterbos fc Kennicutt]ll98Sl ). 
and adopt parameters for the thickness of the disc identi- 
cal to those used in 105, with a constant thickness for the 
disc of 350 pc (which i s roughly the thickness of the MW 
disc, llvezic et al.l J2008j)) out to 16 kpc, at which point we 
assume the stellar disc begins to flare with a scale height 
that increases linearly with radius. We set a maximum scale 
height of 1.69 kpc at radius of 30.5 kpc, beyond which we 
assume that the disc has constant thickness. We then inte- 
grate along the line of sight through this flaring exponential 
disc and project the velocities of objects on circular orbits 
about M31 onto the line of sight. This produces an average 
velocity map for the disc of M31, which we display in Fig. [4] 
Once our disc model has been constructed, we need to 
select a sample of fields from our DEIMOS survey that will 
provide the most reliable kinematic comparison with respect 
to the velocity of the disc. As the disc of M31 is not infinites- 
imally thin, but possesses some unknown scale height, any 
line of sight taken through the galaxy traverses a significant 
depth. Given the inclination of M31, some lines of sight will 
traverse larger depths than others, which could have the ef- 
fect of smearing out the velocities of objects with respect to 
the disc model. This is illustrated in Fig. 8 of 105. They find 
that objects along the major axis of M31 are less susceptible 
to this effect than those that are located off the major axis, 
and therefore we limit our initial study to fields along the 
major axis. 

A further complication in field selection arises from 
MW contamination. Our colour selection criteria means 
that we inevitably observe Galactic K dwarf stars within 
our sample, as these lie coincident with M31 RGB stars in 
the CFHT and INT CMDs. Eliminating these stars from 
our sample is straightforward in the South West (SW) re- 
gion of M31, as the disc of Andromeda and the halo of 
the MW occupy distinct positions in heliocentric velocity 
space. Assuming the Besancon model is a good description 
of the foreground populations in the direction of M31, it 
can be shown that the Galactic population peaks at Vhei=- 
61kms _1 , and the contribution of MW K dwarfs to our sam- 
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Table 1. Properties of fields analysed in this work 
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pie at v he i < — lOOkms 1 is very low (|Robin et alj 12004 
105). Given that the average rotational velocity in the SW 
of M31 is less than -300 kms" 1 (105), we are able to cleanly 
separate M31 stars from MW field stars. However, in the 
North East (NE) of M31, the average heliocentric veloc- 
ity of the M31 disc typically ranges between -100 kms -1 
and -200 kms -1 , resulting in a significant overlap between 
Galactic and M31 populations, making it difficult to distin- 
guish between the two. While it is possible to remove some of 
this contamination by examining the strength of the Sodium 

o 

doublet (Nal), located at a rest wavelength of ~ 8190 A, 
this is not a perfect discriminator. One can also eliminate 
some foreground contamination via a comparison of photo - 
metric and spectroscopic metallicities l|Gilbert et ali l2006). 
but given the uncertainties on the individual spectroscopic 
[Fe/H] of our observed stars (discussed above), we still retain 
a significant population of contaminants within our sam- 
ple. There is also contamination in the NE from M31 sub - 
structure (I05 lChapman et al.ll2006l : iRichardson et al.ll2008l ) 
which can be difficult to separate from the M31 disc in the 
NE. For these reasons, we limit our initial study to the SW 
major axis. We hope to analyse the NE population in a 
future paper. These criteria leave us with a sample of 21 
fields along the SW major axis, highlighted in red in Fig. [1] 
The positions and properties of these fields can be found in 
Table [1] Two of these fields (231Dis and 232Dis) were ob- 
served as part of our ultra-deep M31 disc survey (Chapman 
et al. in prep.) and were integrated for 4.5 hours with the 
600 line/mm grating, which allowed us to make more robust 
measurements of individual stellar metallicities than for our 
other fields. 
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Figure 5. Histograms for both heliocentric (top) and disc lag 
(bottom) velocities of all stars within our sample of 21 fields. Re- 
gions expected to be inhabited by thin disc (light blue), thick disc 
(red), halo (green) and MW foreground (grey) are highlighted. 



To isolate a potential thick disc population in our sam- 
ple of 21 fields, we must first define a statistical measure 
of what constitutes a lagging population with respect to the 
thin and extended discs. We must also ensure that our defini- 



tion is able to distinguish between this population and that 
of the meta l-poor M31 halo which, as it is a non-ro tating 
component (|Chapman et al.| [2006; Kal irai et al.ll2006T ). also 
lags behind the disc. In Fig. [5] we display two histograms, 
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Figure 6. Our initial sample of fields, selected for their position along the South West major axis as described in the text, are shown in 
order of increasing (projected) radius. Gaussian fits indicating the thin and (where applicable) thick disc are shown as magenta and blue 
curves respectively. Our selection criteria are overlaid, with the dashed lines representing our 2a cuts and the solid lines representing our 
Gaussian cuts both of which sample roughly the same region of velocity space. 



one with the heliocentric velocities (vhel) of the stars in our 
21 fields, and one with their velocities with respect to the 
disc (viag) and we highlight the regions we expect each of 
these components to inhabit, along with where we expect to 
see contamination from halo K-dwarfs in the MW. 

We do this using two separate methods. The first is 
to fit Gaussians to both a disc component, located on or 
around vi a9 =0 kms -1 , and a broad halo component cen- 
tered on or around vi ag —-300 kms -1 . We then define a thick 
disc population to encapsulate anything that lags the thin 



disc by > 2a of the thin disc peak value and we implement 
a lower cut on this population by requiring the contribu- 
tion from the halo would be <1 star per velocity bin (20 
kms -1 ), thus minimizing the contamination. In fields where 
there is no obvious halo component to fit to, we use a Gaus- 
sian centered on -300 kms -1 with a dispersion of 90 kms -1 
jChapman et alj|2006h . and normalize it with respect to the 
thin disk by assuming that the halo contributes ~10% to 
the total number of stars within the field (a conservative 
estimate, given that the stellar halo contributes << 10% to 
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the total stellar light in disc galaxies). The second method 
is to fit multiple component Gaussians to each of the fields. 
We apply a Gaussian Mixture Modelling (GMM) technique, 
which allows the number of Gaussians to vary freely be- 
tween 1 and 7 components. To discern which model best 
fits the data, we apply a likelihood ratio test (LRT) to the 
resulting probabilities of the fits. The us e of the LRT in as- 
tronomy was popularised bv lCashl (| 19791 ). and is often used 
in the literature to determine whether the properties of an 
observed stellar population can be wel l described by single 
vs. multiple Gaussi an components fe.g. lAshman et aL l l 1994 
ICarollo et aT1l2O10D . The LRT compares the likelihoods of 
nested models (in our case, a mixture of Gaussian compo- 
nents) to determine whether applying a model with addi- 
tional parameters produces a significantly better fit than a 
simpler model. This is done by calculating the LRT statistic, 
— 2\n(£i / C2) , where C\ and £2 represent the likelihoods of 
the simple and complex model respectively, and comparing 
it with a x 2 distribution with degrees of freedom equal to 
the difference between the number of parameters in the two 
models (3 in our case). For a model with additional param- 
eters to be accepted as a statistically better fit, this ratio 
must be greater than 7.82 which corresponds to a P-value 
of < 0.05. In general, this technique converges on fits with 
three components (a thin disc, a halo and a thick popula- 
tion) though there are a few exceptions. We shall discuss 
these fits in greater detail in the following section. Where 
this technique converges on fits that identify a lagging com- 
ponent that is distinct from both the thin disc and halo, we 
define a sample of highly probable thick disc stars by ap- 
plying a standard Bayesian classification scheme to assign 
each star a probability of being a member of the thin disc 
P(thin), thick disc, P(thick), or halo, P(halo), population 
based on their velocity, and the properties of the Gaussian 
fits to each population on a field by field basis. We define a 
star as being a highly probable member of the thick disc if 
P(thick)> 0.997. The results of both these techniques can be 
seen in Fig. [6] The velocity cuts for stars selected using our 
2a technique are shown as dashed lines, and the range of ve- 
locities selected using the Bayesian classification technique 
are marked with solid lines. It can be seen that both tech- 
niques isolate a very similar population. In Fig. we plot a 
CMD showing the V-I colours of the thin (blue points) and 
thick (magenta po ints) populations, and we overlay Dart- 
mouth isochrones (|Dotter et al.ll2008T ) of [cv/Fe] = +0.2 and 
an age of 8 Gyrs (in line wi th the estimated r ange of ages 
for the thin disc of 4-8 Gvrs. lBrown et aill2006l ) with metal- 
licities ranging from [Fe/H] = -0.4 to [Fe/H] = -1.5. Both 
thin and thick populations inhabit roughly the same region 
in this CMD, which we shall discuss in more detail in §4.3. 

Finally, we note that in both selection methods, we ex- 
pect some cross contamination between the thin and thick 
disc components, as the two populations significantly over- 
lap. However, we assume this contamination will be lower 
in our cuts based on the Gaussian fits, as these are more 
conservative. Therefore, we use these cuts predominantly 
in this paper when referring to clean thin and thick disc 
samples. We have also assumed both components have sym- 
metric, Gaussian distributions in velocity, which may not be 
the case, and this could cause further contamination if the 
populations are skewed. We also expect some contamination 
from the halo, however, given that the disc is the dominant 




(v-i) 



Figure 7. CMD for our thin (blue points) and thick (magenta 
points) populations in standard Landolt V and I colours. Dart- 
mouth isochrones with [a/Fe] = +0.2 and an age of 8 Gyrs are 
overlaid with metallicities from left to right of [Fe/H]=-1.5, -1.0, 
-0.5, -0.4. The colours of the remaining stars in our DEIMOS 
survey are also plotted in light blue. 



population in all our fields, we expect this contamination 
to be negligible in comparison to the cross contamination 
between the discs. 

3.2 Testing the statistical significance of our 
sample 

Before we analyse our sample, we test the significance of 
our thick population to ensure it is not merely consistent 
with noise above a thin disc population plus smooth halo 
component. To do this, we fit a single Gaussian to the disc 
and halo components, as described above, then calculate the 
deviation of the data from the fit for all velocities greater 
than the peak disc velocity (i.e. the right hand side of the 
disc fit), normalizing it to the expected contribution from 
the Gaussian in this region. We then define the noise to be 
1.5 times the median absolute deviation of this sample. We 
repeat this exercise for all velocities less than the thin disc 
peak and greater than -200 kms -1 in the lag frame, in this 
case comparing to the expected contribution from both thin 
disc and halo fits. This allows us to work out the signifi- 
cance of our thick disc population, a con f. In all cases where 
the GMM identified a thick disc component, we find that 
our excess above the thin disc plus halo model has a signif- 
icance of > 3a (see table [2}. For the fields where the GMM 
converged on a 2 Gaussian fit (232Dis, 222Dis, 107Ext and 
wllold), we find a con f < 3a. We also identify an additional 
3 fields (lOlDis, 166Dis and 227Dis), where a con f < 3a. 
Two of these fields are located at radii of ~ 15 kpc, where 
there may be residual contamination from the bulge compo- 
nent. This may also explain the large dispersions (of order 
50kms _1 ) seen in our innermost fields. Excluding these, we 
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Figure 8. Top panel: The difference in velocity, Av, between 
the thin disc and thick component as a function of projected 
radius. This lag appears to be approximately constant as a func- 
tion of radius, with an average lag of 46.0 ± S^kms" 1 . There 
appears to be a slight in crease in lag in the outermost part, how- 
ever this is largely driven by fields that lie off the major axis of 
M31, and therefore the velocities are less reliable. Middle panel: 
Dispersion, a v , of both thin (black squares, solid line) and thick 
(red triangles, dot-dashed line) components are plotted as a func- 
tion of projected radius. The thin disc appears to maintain a 
constant dispersion of a t f lin =35.7H.O kms^ 1 , however the thick 
component appears to decrease somewhat at larger radii. Bot- 
tom panel: Average spectroscopic metallicity of thin and thick 
components as a function of projected radius. Neither component 
evolves with radius. 



are left with 14 of our 21 (2/3) fields where we confidently 
detect a thick component. We shall focus on these fields in 
the remainder of our analysis, but we shall discuss the sig- 
nificance of the non-detections in 85. 



4 RESULTS 

4.1 Kinematic and structural properties of the 
thin and thick discs 

In this section we present measurements for the kinematic 
and structural properties of the thin and thick discs. Prop- 
erties of individual fields can be found in Table [2] while the 
average properties for both components can be found in Ta- 
ble El 




Figure 9. Top panel: Plot of the density of stars (N*/arcmin 2 ) 
in the thin disc against R pro j. The densities are calculated by 
first separating our sample by their target prioritisation (A or B, 
see §3), then counting all stars with v thin > v lag > v thin + 2a thin 
and multiplying these values by 2 (i.e. assuming the distributions 
are symmetric) for both prioritisations. We then calculate p* == 
n s nt/n — n^, and combine these results from priority A and B. 
Fitting an exponential profile to these points we deduce h r = 
7.3 ± 1.1 kpc. Solid line represents the best fit to the data from 
a weighted-least-squares routine, and the shaded region indicates 
the lcr errors from the fit. Bottom panel: As above, for the thick 
disc. Here we count all stars with v thick - 2a thick > v tag > v thick 
and multiply by two again. Fitting an exponential profile to these 
points we deduce h r = 8.0 ± 1.2 kpc. 



4-1.1 Velocity lag and dispersion profiles 

In this first section, we initially address the thin and ex- 
tended discs of M31. In 105, the extended disc was identified 
as a stellar disc that, while appearing in many respects to 
be similar to the classical thin disc, was a separate entity 
that was clumpy in terms of its structure, and lagged be- 
hind the classical disc in terms of its kinematics. As we are 



limiting our study to one slice down the major axis of M31, 
we do not attempt to comment on the global 'dumpiness' 
of this extended disc, but we return to the issue of the ve- 
locity lag and distinction from the thin stellar disc. As we 
have analysed the disc frame velocities for all our fields us- 
ing a rotation curve that differs from the one used in 105, 
it is useful for us to determine whether the increasing lag 
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Figure 10. Here we compare the scale lengths of the thin and thick discs in M31 to those of other galaxies with observed thick discs in 
order to infer their scale heights. Left panel: Here we plot thin disc scale lengths against thick disc scale lengths for 34 external galaxies 
(YD06) as black points, and overplot the same measurements for the MW (blue triangle). We fit a linear relation to these points with a 
gradient of 1.3. We plot our result for M31 as a red square, and it is in excellent agreement with this relation. Centre panel: Here we 
plot scale height, zo, against scale height for the thin disc of the YD06 sample plus the MW and derive a best-fit linear function with a 
gradient of 0.17. From this, we can estimate a scale height for the M31 thin disc of 1.1±0.2 kpc, which we overplot on the relation as a 
red square. Right panel: We now plot the same, but for the thick disc and find a best fit gradient of 0.35, and therefore infer a scale 
height for the thick disc in M31 of 2.8±0.6 kpc, which is overplotted as a red square. 



with respect to the classical thin disc is seen here also. In 
105, they split their sample of 21 fields into an inner (with 
Rproj < 20 kpc) and outer (with 20<R pro j < 30 kpc) sam- 
ple to determine the average properties of the disc and ex- 
tended disc. For their inner (classical) disc sample, they cal- 
culated an average velocity for the disc in the disc lag frame 
of viag = — 17.0 kms -1 and a dispersion ofcr„ = 50.0 kms . 
In the outer (extended) sample they calculated an aver- 
age velocity of vi ag — — 16.0 kms -1 and a dispersion of 
a v = 51.0 kms -1 . If we perform the same analysis for our 
study, we find an average lag of vi ag = — 14.8 kms -1 for our 
inner fields and vi ag — —25.5 kms -1 for our outer fields. 
However, we note that this value is calculated with the in- 
clusion of fields 107Ext, wllold and 167Hal, which have 
very large lags of vi ag < — 55 kms -1 compared with the 
other fields. We note that these fields are located slightly 
off the semi major axis (see Fig. [T]) where our interpo- 
lated disc-frame velocities are subject to larger uncertain- 
ties. If we exclude these fields, we find an average lag of 
viag = — 14.9 kms -1 , very similar to our inner sample. We 
therefore conclude that there is a negligible difference in the 
lags of the classical and extended disc behind circular veloci- 
ties. For these samples we also calculate average dispersions 
of a v — 42.7 kms -1 and a v = 30.0 kms -1 , implying that 
the extended disc has a lower dispersion than the classical 
disc. However, in our inner sample, we are more likely to 
see residual contamination from the bulge and we also have 
a large proportion of fields for which we could not cleanly 
isolate the thick disc (~ 40% cf. ~ 20% in the outer sample). 
These factors may cause us to overestimate the dispersion 
of the disc in these regions. From these results, we there- 
fore find no concrete reason to assume that the extended 
disc is a separate component from the classical disc and we 
treat these two components as one thin stellar disc in the 
remainder of our analysis. 

By using the information from our Gaussian fits to 
the thin and thick components, we can comment on their 
global kinematic properties, and discuss any variation of 



these properties with radius. In Table [5] we show the peak 
velocities and velocity dispersions of both thin and thick 
(where applicable) components in each field, with associ- 
ated errors from the GMM fits. Where both thin and thick 
components are detected, we compute the lag between the 
two components, Av — Vthin — v t hick , and plot this lag as a 
function of radius in Fig. [5] The 14 fields for which a thick 
disc component is reliably detected cover a range of radii 
from 15.2 to 39.6 kpc. In the top panel of Fig.fS] we can see 
that the lag between the two components does not appear 
to increase with distance from the centre of M31, and shows 
an average lag of (Av) — 46.0 kms - . 

We also plot the dependence of velocity dispersion, athin 
and (Jthick, for both components with radius in the mid- 
dle panel of Fig. [5] For the thin disc, we fit both a con- 
stant relation and a single power law to the data. The lin- 
ear power law suggests a decrease in dispersion with ra- 
dius, with a gradient of -0.87kms -1 kpc -1 , however this 
fit is not statistically better than a constant fit, with an 
average dispersion of athin = 31.6 ± 1.1 kms -1 (reduced 
X 2 of 5.6 vs. 5.2). We do the same for our thick disc re- 
sults, and we find that a linearly decreasing profile where 
vthick = — 0.8(±0.2)i? pro j + 66.1(±5.8) has a marginally 
better fit to the data than a fit with no evolution, how- 
ever the difference in negligible (reduced \ 2 °f 1-2 vs. 1.4), 
and deemed insignificant in a \ 2 significance test. Even if 
we were to accept this fit as preferred, we note that the two 
outermost fields situated at 34.2 and 39.6 kpc, are perhaps 
the driving force in the decreasing dispersion seen in our 
thick disc component. As these field are the furthest out 
in our survey, they also suffers from the greatest chance of 
halo contamination in our sample, and therefore could be 
unreliable. If we exclude these final points from the fit, we 
find that athick is best fit with no evolution as a function of 
radius, with an average dispersion of 50.8±1.9 kms -1 . We 
therefore conclude that our data cannot tell us anything re- 
liable about the dependence of these kinematic properties 
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Figure 11. Here we display the spectroscopic MDF for all stars in our thin and thick components (shown as blue hatched and red 
solid histograms respectively) In the left panel, we show the MDF using all stars for which metallicities can be reliably measured (i.e. 

o o o 

S:N> 10 A -1 , and the middle and right panels apply lower quality cuts (S:N> 8 A -1 and S:N> 5 A -1 ). For our lower S:N cuts, we note 
that the median [Fe/H] values for both populations remain similar, and the dispersions (inter-quartile range) begin to increase, losing 
some of the detail of the shape of the MDF. In all cases the median [Fe/H] of the thick disc is more metal poor than the thin by ~ 0.1 
dcx. 
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Figure 12. Composite spectra of both the thin disc and thick 
component, using both 2cr (top panel) and Gaussian (lower panel) 
cuts to isolate the thick component. The composite spectra are 
constructed from stars within the selection regions that possess 

o 

a S:N > 3.0 A -1 . We find that the average metallicity for the 
thin component is more metal rich than the thick by 0.2—0.3 dex. 
These results are consistent with composites formed from spectra 

o 

with S:N < 3 A" 1 and with our ficld-by-ficld metallicity estimates 
(Fig. [§}. We also show the locations of a number of Fe I lines 
present in these spectra. 



with radius, and allow us to merely calculate the average 
kinematics of both components. 



4-1-2 Scale length of the thin and thick disc 

To determine the scale lengths of our two disc components, 
we need to calculate the number density of thin and thick 
disc stars within our DEIMOS field of view. There are 2 
complications we must consider before we proceed. Firstly, 
the two components are not completely distinct from one 
another, and in all fields, we observe some overlap. Secondly, 
owing to our selection criteria (discussed in § 3), we prioritise 
stars of certain colours and magnitudes above others, and 
this must be considered when calculating densities on a field- 
by-field basis. 

We determine the number of stars associated with the 
extended thin disc, n s , in each of our fields by integrating 
the Gaussian we have fit to this component. To determine 
the density of stars contained in the thin disc, we multiply 
n s by the total number of available target stars within our 
DEfMOS field that fall within our selection criteria, n t , and 
divide this by the total number of stars that were observed 
with our DEIMOS mask, n . We then subtract the density 
of background stars rib, which is computed from a number of 
fields on the edge of our survey region; i.e. p* = n s n t /n — rib- 
To account for our prioritised selection technique, we per- 
form this calculation separately for our priority A and pri- 
ority B stars, then combine these measurements. We repeat 
this calculation for the thick disc. We plot the results in 
Fig. [9] where we apply a weighted least-squares exponential 
fit to our data points, and determine h r — 7.3 ±1.1 kpc for 
the thin disc and h r = 8.0 ± 1.2 kpc for the thick. Com- 
paring this to previous calculations for the scale length of 
the thin and extended discs, we find that the extended disc 
has a larger scale length than the exponential thin disc, 
(5.1 ± 0.1 kpc, 105). The value of 7.3 kpc that we derive 
is slightly higher than that derived in 105 of 6.6±0.4 kpc, 
and with much larger error bars but the two are consistent 
within their la uncertainties. The difference between the 
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two values can be attributed to two factors. Firstly, in 105, 
they included fields from the NE of the galaxy, plus fields 
located away from the major axis, where we have have sam- 
pled fields solely from the SW major axis. Secondly, in 105 
they did not fully address any biases that may have been 
introduced by our two-tiered prioritisation system. Finally, 
we note that the thick disc appears to be more radially ex- 
tended than either the thin or extended disc, although it is 
consistent with the scale length of the extended disc within 

its l<T-errors. 

In previous work lYoachim fc Dalcantonl (120061 ) (here- 
after YD06) measured the scale lengths of 34 edge-on disc 
galaxies using a photometric fitting technique, and found 
that the scale lengths of the thick discs were larger than 
those of the thin discs by a factor of~ 1.3. We plot their re- 
sults in the left panel of Fig. 1101 and overlay a linear relation 
with a gradient of 1.3. We add to this our results for M31, 
using an average value for the thin disc from the range of 
scale lengths derived for the thin and extended discs (5.1- 
7.3 kpc) of 6.3 kpc, and our calculated val ue of 8.0 kpc. We 
also overplot the result for the MW (using [juric et al1l2008l 
values of 2.6 and 3.6 kpc for thin and thick discs respec- 
tively), and note that M31 sits in excellent agreement with 
this relation. 



4-1.3 Inferring the scale heights of the thin and thick discs 

Owing to the inclination of M31, we are unable to mea- 
sure the height of either the thin or thick disc components 
directly. No photometric excess above a typical bulge or ex- 
tended disc profile is observed when performing minor axis 
star counts (llrwin et al.ll2005T) . suggesting that these com- 
ponents dominates the surface profile out to large radii. In 
order to infer probable scale heights for both components, we 
make use of the properties of the 34 edge-on galaxies mea- 
sured by YD06. As the scale lengths and heights of both 
thin and thick discs in each of these galaxies were derived, 
it is possible for us to search for a relation between the scale 
length, h r , and scale height, zo of each component. In the 
central panel of Fig. 1101 we pl ot h r vs. zp for t he YD06 
sample as well as for the MW (jlvezic et al.l 1200*81 ). and fit 
it with a linear relation, on which we force an intercept of 
(0,0). We find that the data are well fit with a gradient for 
this relation of 0.18±0.04, though there is significant scatter 
beyond ~ 9 kpc. From this, we deduce zo=l-l±0.2 kpc for 
the M31 extended disc (using h r =7.3 kpc). We repeat this 
for the thick disc (shown in the right panel of Fig. 1 10 fl and 
find that these values are well fit with a linear relation of 
gradient 0.35±0.06, giving us z =2.8±0.6 kpc for the M31 
thick disc. If these values are correct, then not only are the 
discs of M31 more radially extended than those of the MW 
by a factor of ~ 2 — 3, they are also significantly thicker. 



4-1-4 Contrast of the thin and thick discs 

In the previous sections, we have derived the density in each 
field of both our components as a means to determine the 
scale lengths. We now use these densities to work out how 
much of the total (disc related) stellar population is con- 
tained within either component. There are several caveats 
to such a comparison that should be mentioned. Firstly, our 



sampling of the field is likely to have an effect on our field- 
to-field estimates of the stellar density (which we discuss 
further in § 5). Secondly, as the disc is not observed edge- 
on, we are measuring a 2D projection of the densities which 
is difficult to interpret. We also note that the measurement 
errors associated with the densities of each field (shown in 
table [4} are significant (of order ~ 50%). 

With this in mind, we find that, on average, the thick 
disc component accounts for 35% of the total stellar density, 
with an inter-quartile range of ±10%. In the Milky Way, we 
know that the thick disc contributes to ~ 10% of the stellar 
density in the solar n eighbourhood, and accounts for ^1/3 
of the total disc mass (I Juric et al.| [2008; Sch onrich fc Binnevl 
l2009al ). comparable to what we derive here. 

From our calculated contrasts and individual density 
profiles for the thin and thick discs, we can estimate the 
mass contained within the thick disc component using val- 
ues for the mass of the thin disc from the literature. From 
our analysis above, we have determined that the thick disc 
contributes 35±10% of the total stellar density, meaning 
the thick:thin disc density ratio is of order 55±15%. We 
can also estimate this fraction by integrating our stellar 
density profiles (Fig. [9} over the limits of our data, and 
from this we calculate a thick:thin disc density ratio of 
~65%, which is in good agreement with our contrast es- 
timate. If we assume that both discs are composed of sim- 
ilar stellar populations, we can set the mass rati o between 
the di sc to be equivalent to the density ratio. In I Yin et al.l 
(2009), they quote a total mass for the thin stellar disc of 
M*,thin = 5.9 x 10 10 M q, c alculated from the m ass models 
of lWidrow et al.1 (120031 ) and lGeehan et al.1 (|2006l ). From this 
we estimate that the total mass of the M31 thick disc lies in 
the range 2.4xl0 10 M Q < M t>thick < 4.1 x 1O 1O M . As we 
are unable to determine the full radial and underlying lumi- 
nosity profile for the thick disc, these values are obviously 
prone to large errors introduced by our simplifying assump- 
tions, and the mass of the thick disc may be lower than our 
quoted range. For example, if we just compare the number of 
stars we detect in the thin disc throughout our entire sample 
with the number we detect in the thick disc by integrating 
the fitted Gaussians in Fig. [6] we find a thick:thin disc ratio 
of 20%. If we assume this value for the ratio of the masses 
between the components, our lower limit on the thick disc 
is reduced to M* jt h in = 1.2 x 10 10 Mq. A future study of 
the thick disc which includes fields from the entirety of our 
survey will help us to better constrain both the radial profile 
and mass of this component. 



4.2 Spectroscopic metallicities 

In this section we present the spectroscopic values of [Fe/H], 
both for individual stars, and for the composite spectrum of 
each component. Measuring individual metallicities from the 
Ca II triplet for the stars in our survey, with S:N typically 

o i 1 

betwe en 5-15 A -1 , is quite problematic. In lBattaglia et al.l 
(2008), they show that the 'best case' errors in measuring 
the equivalent widths of the Ca II lines, AEW, scale with 
S:N as: 



AEW = 



Vl.5 x FWHM 
S : N 



(3) 
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Table 2. Kinematic properties of thin and thick disc components 



Field 


Vthin,lag ( km S 1 ) 




Vthick,lag (kms 1 ) 


o- tMcfc (kms x ) 


^conf 


[Fe/H] tMn 


{Fe/H] t hick 


^ZiOJ— 'la 


6.8±5.0 


c:c 9-LQ 9 

OO . £^ZO.tZ 


-141,9+7,5 


41.2+11.8 


37.0 


-0.8+0.1 


-1.0+0.1 




q 7_|_ in c 


uo. / zno.o 


N /A 


N /A 


9 8 
z .0 


7-(-n 9 
-u. / znu.z 


N /A 


166Dis 


-7.5+6.9 


48.9±4.6 


N/A 


N/A 


0.5 


-0.8+0.2 


N/A 


1 OfiFlic; 


1 1 fi-LQ 9 




qc; n_l_Q 9 


^9 ^A-A 1 


0.0 


n Q-i-o ^ 
-u.yznu.o 


-r .uinu.o 




1 R AA-7 9 


Q9 fl-t-A 1 


^4 8-1-8 7 


c:i n-t- ^ 9 
or .uzno.z 


1 n 


7-1-0 9 
-u. / znu.z 


1 0-1-0 9 


994T)iq 

ZiZi'i.L/lo 


-34.4+6.2 


94 q-|-9 c; 




r:c: n_|_l 9 r: 


5 1 


-0.9+0.1 


-1.1+0.2 


232DiS 


-37 2+3 




N/A 


N/A 


2.0 


-0.8+0.1 


N/A 


104Dis 


o ( . y zlu . <j 






CO QJ-7 9 


9 


-0.7+0.2 


-0.9+0.1 


220Dis 


-12.5±5.0 


39.6±2.9 


-64.8+12.1 


54.7+10.9 


26.7 


-1.1+0.3 


-1.6+0.3 


213Dis 


-10.9±3.2 


31.9±6.3 


-80.5+10.0 


55.0+8.8 


39.2 


-0.8+0.3 


-1.5+0.3 


102Dis 


-18.8±7.1 


22.2±5.2 


-61.5+11.6 


48.1+5.6 


47.1 


-0.9+0.2 


-1.0+0.2 


231Dis 


4.7+15.2 


22.9±4.2 


-35.9+7.8 


53.4+7.2 


43.4 


-0.9+0.1 


-1.1+0.2 


223Dis 


-8.1+5.9 


20.1±4.9 


-55.7+9.9 


51.1+6.8 


14.7 


-0.7+0.2 


-1.0+0.3 


lOlDis 


-16.2±3.9 


35.8±5.8 


-57.3+8.6 


44.9+6.2 


3.1 


-1.1+0.2 


-1.3+0.3 


222Dis 


-25.9±4.8 


33.6±2.9 


N/A 


N/A 


0.8 


-0.7+0.1 


N/A 


221Dis 


-27.5±5.5 


35.0±3.6 


-121.2+7.4 


52.2+6.3 


27.1 


-1.1+0.2 


-1.4+0.2 


50Disk 


-16.9±8.1 


34.2±5.9 


-149.2+13.2 


42.3+7.2 


12.6 


-1.0+0.2 


-1.2+0.2 


107Ext 


-55.2+6.9 


24.8±4.4 


N/A 


N/A 


2.4 


-1.0+0.3 


N/A 


wllold 


-55.1+10.2 


28.8+5.6 


N/A 


N/A 


1.4 


-1.0+0.3 


N/A 


167Hal 


-72.6±8.3 


22.9±5.4 


-120.3+12.2 


35.1+7.8 


23.5 


-0.9+0.3 


-1.0+0.3 


148Ext 


-17.2±7.0 


41.3±6.6 


-154.5+11.0 


25.7+5.1 


28.2 


-0.9+0.3 


-1.0+0.2 



Table 3. Average properties of thin and thick disc components 
derived in this work 



Component <r„(kms 1 ) h r (kpc) zo (kpc) [Fe/H] S j, ec 

Thin disc 35.7+1.0 7.3+1.1 1.1+0.2 -0.7+0.05 

Thick disc 50.8+1.9 8.0+1.2 2.8+0.6 -1.0+0.1 



Table 4. Field by field densities of thin and thick disc stars 



Field 


Ptthin (*/arcmin) 


ptthick (*/arcmin) 


228Dis 


57.6+22.3 


N/A 


227Dis 


56.9+22.3 


N/A 


166Dis 


20.4+9.4 


N/A 


106Dis 


36.2+15.3 


23.6+10.5 


105Dis 


45.3+18.6 


41.0+16.3 


224Dis 


20.1+8.4 


17.1 ±7.1 


232DiS 


35.9+15.7 


N/A 


104Dis 


35.7+14.6 


6.6+4.5 


220Dis 


19.9+8.9 


13.5+6.4 


213Dis 


23.0+12.2 


9.8+3.0 


102Dis 


24.0+12.8 


8.9+4.5 


231Dis 


13.0+6.6 


16.6+8.1 


223Dis 


15.1+6.7 


12.2+5.6 


lOlDis 


26.5+12.9 


11.1+4.8 


222Dis 


4.9+ 2.1 


N/A 


221Dis 


14.4+6.4 


4.8+3.3 


50Disk 


8.1+3.9 


2.0+0.6 


107Ext 


9.6+7.8 


N/A 


wllold 


9.0+6.5 


N/A 


167Hal 


1.3+0.8 


0.8+0.5 


148Ext 


1.1+0.9 


0.7+0.5 



assuming no contamination from residual sky lines and no 
covariance noise, where FWHM is the full width at half- 

o 

maximum of the CaT lines which is typically 2-3 A. Using 
this equation, we can determine the average errors in [Fe/H] 
for stars in our sample at different S:N, and we find that 

o 

for spectra with S:N of (5, 8, 10) A , the errors in their 
calculated metallicity are of order (0.6, 0.4, 0.3) dex. To 
demonstrate the effects of these large errors on the metallic- 
ity distribution function (MDF) of our sample, we present 
histograms of the individual spectroscopic metallicities for 

o 

three quality cuts, one at S:N> 10 A -1 (left panel), one with 

o o 

S:N> 8 A" 1 (middle panel) and one with S:N> 5 A" 1 (right 
panel) in Fig. 1111 In both plots, the blue hatched histogram 
represents our thin disc sample and the filled red histogram 
represents our thick disc population. For our higher qual- 
ity spectra, we calculate a median metallicity for the thin 
disc of [Fe/H]=-0.9 with a dispersion of 0.4 dex (from the 
inter quartile range, IQR). For our thick disc population we 
calculate a median of [Fe/H]=-1.0 also with a dispersion of 
0.4 dex. We note that the distributions of both populations 
deviate from a Gaussian distribution, with a kurtosis of - 
0.6 and -0.3 for thin and thick discs respectively, implying 
a broad peaked distribution, with narrow tails. Both distri- 
butions are skewed towards lower metallicity with skewness 
a — —0.4 and a = —0.3 for thin and thick disc. For our low- 
est S:N cut, however, much of this information is lost. While 
the median [Fe/H] remains very similar with [Fe/H]=-1.0 for 
the thin disc and [Fe/H] =-1.1 for the thick, the distributions 
begin to broaden, with dispersions of 0.5 dex for both popu- 
lations, and present almost no skew (a = —0.1 and a = —0.2 
for thin and thick disc). This shows that by including data 
with larger measurement errors, we wash out our MDF con- 
siderably, and lose any meaningful information. As a sanity 

o 

check, we compare the MDF for all stars with S:N> 10 A -1 
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with one for stars with S:N> 15 A -1 , and find both the me- 
dian values of [Fe/H] and general distributions to be com- 
parable. We note that by requiring such a high S:N cut on 
our individual measurements of [Fe/H], we bias our sample 
towards more metal-rich stars as these will have intrinsically 
stronger Ca II lines. 

Owing to the large errors associated with these mea- 
surements, this analysis provides a crude indication of the 
metallicities of both discs, and so to get a more accurate 
estimate of the average metallicities of both populations 
we construct composite spectra for both components (using 
both the 2a and Gaussian velocity cuts) by co-adding the 

o 

individual spectra of all stars with S:N> 3 A -1 , weighted 
by their S:N values. The resulting S:N of the composite is 
much greater than the individual spectra (S:N~ 60— 100 cf. 
S:N~ 3 — 25), allowing a better fit to the Call lines. We use a 

o 

cut of S:N> 3 A -1 as below this the velocity uncertainties of 
our stars begin to significantly increase (as discussed in § 3). 
As we shift all spectra to the rest frame before co-adding, in- 
cluding spectra where the velocity is uncertain could smear 
out the Ca II lines, resulting in an over-estimate of [Fe/H] 
for the composite. We note that the results from our com- 
posites are only indicative of an average metallicity for each 
component, and can tell us nothing about the metallicity 
dispersion for the discs. We display the resulting composites 
in Fig. 1121 The top 2 panels show the thin and thick spectra 
for the 2a velocity cuts, while the bottom 2 panels show the 
same, but for the Gaussian velocity cuts. In the case of the 
2a cuts, our thin composite comprises 511 stars that match 
our kinematic and quality criteria, while our thick composite 
is constructed from 78 stars. For our Gaussian cuts, these 
numbers fall to 380 and 52 stars respectively. We find an off- 
set of order 0.2 dex between the thick and thin components 
for the 2a cuts, with the thick disc being more metal poor at 
[Fe/H]= -0.9 ±0.1 compared with [Fe/H]= -0.7 ±0.05 for 
the thin, inconsistent within their respective la errors. For 
our Gaussian cuts, we find the thick disc to be more metal 
poor, giving us a larger difference in metallicity between the 
two components of 0.3 dex (with [Fe/H]= —1.0 ±0.1 for the 
thick disc compared with [Fe/H]= — 0.7±0.05 for the thin), 
although the two results for the thick disc are consistent 
within their la errors. We also note that we are liable to ex- 
perience non-negligible thin disc contamination of our thick 
disc component, which could cause us to over estimate the 
average [Fe/H], so the true difference could be larger still. 
We note that these results are consistent with performing 
the same analysis on composites constructed from spectra 

with S:N> 10 A" 1 . 

Finally, the continuum fit to the third line in our com- 
posite spectra, particularly for our thick disc selection, gives 
us some cause for concern. Could this metallicity difference 
we derive be driven by poor continuum fitting in this region 
of the spectrum? To investigate this, we analyse the [Fe/H] 
for the thin and thick discs again, using solely the first two 
lines (Call8498 and Call8542- In the case of our simple 2a 
cut, this narrows our difference in metallicity slightly from 
0.2 dex to 0.15 dex, with [Fe/H]= -0.85 ± 0.1 compared 
with [Fe/H]= -0.7 ± 0.05 for the thick and thin discs re- 
spectively. However, in the case of our Gaussian cuts, which 
are arguably less affected by cross contamination between 



the components, the metallicity difference of 0.3 dex per- 
sists. 

We also perform this composite analysis on a field-by- 
field basis. The results of this analysis, shown in Table [2] are 
again, less accurate than our overall composite, but they 
suggest a similar offset in metallicity exists in the thin and 
thick components in each field. We plot this result as a func- 
tion of radius in the lower panel of Fig. [8] We find no evi- 
dence for any evolution of metallicity with radius. 

A slight concern in ascertaining the metallicity of a pop- 
ulation from a composite spectrum arises from inaccuracies 
in the estimate that come from combining spectra with dif- 
ferent effective temperatures and l/-band magnitudes, as the 
derived metallicities are weakly dependent on the apparent 
F-band colours of the stars. The rms dispersion in the V- 
band magnitudes within our sample are small (< 0.5 mag 
for both thin and thick discs) as we are sampling only a 
small region of the tip of the RGB, so the error introduced 
by this effect will be very small. However, to further assess 
this, we separate our thin and thick disc spectra into bins 
of 0.2 mags in the V-band and create composite spectra for 
each bin, measuring the metallicity of each. We show a sam- 
ple of these spectra in Fig. 1131 labelled with the metallicity 
and average V-band magnitude. The typical errors in metal- 
licities determined for these composites ranges from 0.1-0.3 
dex. What we see is that the composite thick disc spectrum 
in each bin is more metal-poor than the corresponding thin 
disc composite. We also find that the average metallicities 
for both thin and thick discs agree with those that we de- 
rived from the composites for the entire sample. 

4.3 Photometric Metallicities 

We inspected the photometric metalli cities of our sampl e 
using the Dartmouth isochrone models (iDotter et al. 2008). 
We se lect an age of 8 Gyrs as the work of iBrown et al.l 
(2006) suggests that the age of the disc in these outer re- 
gions varies between 4 and 8 Gyrs. We use an a-abundance 
of [ a/Fe]=±0.2 as it has been shown in va rious works 
fe.g. lReddv et al.ll2006l ; lAlves-Brito et alJl201Ch that the a- 
enhancement of thin or extended stellar disc populations 
typically ranges between [a/Fe]=±0.0 and [a/Fe]=±0.2. We 
then interpolate between these isochrone models for every 
star within our sample to determine its metallicity. We can 
then compare the MDFs for our thin and thick disc sample, 
selected by both the 2a and Gaussian cuts discussed above. 
The results of this are shown in the left panel of Fig. 1141 
This figure shows us that when using this set of isochrones, 
the MDFs of both populations trace each other remarkably 
well. We calculate a median metallicity for each component 
and find [Fe/H] tMn = -0.79 and [Fe/B] thick = -0.80, both 
with IQRs of 0.2 dex. Neither population has a Gaussian 
distribution, with positive kurtosis of ±2.2 for both MDFs 
(i.e. more peaked, with broader tails), and both populations 
are skewed towards lower [Fe/H] with a ~ —1.2 for both 
discs. From this analysis, one might conclude that the two 
discs are chemically indistinguishable. This is in contrast to 
our findings from the combined spectra in §4.2 where we find 
an offset in the average metallicities of thin and thick com- 
ponents of 0.2 dex. As our photometric data are not deep 
enough to detect the MSTO of these fields, we are exposed to 
the age-metallicity-[a/Fe] degeneracy problem. If we analyse 
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Figure 13. Composite spectra of our thin (left panel) and thick (right panel) disc samples, binned in V-band magnitude. Each bin spans 
0.2 mags. We see that in each case, the thick disc is more metal-poor than the thin disc by ~ 0.2 dex. The errors in the values of [Fe/H] 
for these composites ranges from 0.1-0.3 dex. 



our data with isochrones of different ages and abundances, 
we find that the individual metallicities we measure change. 
Increasing the age by 2 Gyrs has the effect of decreasing 
[Fe/H] of a star by ~ 0.05 dex on average (shown in the 
centre panel of Fig I14p and increasing the abundance from 
[a/Fe]=+0.2 to [a/Fe]=+0.4 reduces [Fe/H] by ~ 0.1 dex, 
(right hand panel, Fig. 1 1 4 p . The dispersions, kurtosis and 
skew remain largely unchanged by these variations. These 
findings demonstrate that it may be difficult to discern slight 
differences in metallicity (such as the 0.2 dex measured in 
§4.2) using photometric isochrones without knowing the ages 
and/or a-abundances of the thick and thin disc. Studies of 
the thin and thick discs in the MW have shown that the 
thick disc is both older and more Q-enriched th an the thin 
disc l|Reddv et al.ll200d ; lArves-Brito et alj|2010h . and many 
of the formation scenarios of thick discs suggest this could 
be true for thick discs in general, including M31. Such dif- 
ferences would certainly affect our derived values of [Fe/H] 
for both discs. 



5 DISCUSSION 

In this section we discuss our findings, and comment on 
the morphology of this thick component. First we compare 
our findings with an expected thin+thick disc population 



inclined to us along the line of sight by 77°, by creating a 
model of a galaxy with a thin/extended disc with similar 
properties to those of M31 that has an additional thick disc 
component and analysing it in the same way as our data. 
We then compare the M31 t hick disc to the MW and the 
lYoachim fc Dalcantonl l|2006l ) sample of thick discs. Finally 
we comment on the possible formation mechanisms for this 
component. 



5.1 Comparison with thin + thick disc model 

To lend confidence to our defining the lagging component 
we isolate in the above analysis as a thick disc, we create a 
simple kinematic model of a galaxy with properties similar 
to those of a MW-type galaxy, which has both a thin and 
thick stellar disc, and analyse this in the same manner as 
our data. This is done as follows; first, we create a thin stel- 
lar disc of 9xl0 6 stars, randomly generating radii for each 
assuming the stars are distributed in an exponential disc 
with a scale length equal to that of M31's (6.6 kpc, 105). We 
assign each particle with a velocity randomly drawn from 
a Gaussian population centred on Okms -1 with a velocity 
dispersion of 25kms _1 in the disc frame. We repeat this for 
our thick disc component, assuming a thimthick disc den- 
sity ratio in M31 that is equal to that measured in the solar 
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Figure 14. Photometric MDFs derived from Dartmouth isochrones of varying age and [a/Fe] l lDotter et alj|2008f ) for the thin and 
thick components (shown as blue hatched and red filled histograms respectively) as defined by our 2a cuts. Left panel: Analysis of 
[Fe/H] for thin and disc using [cr/Fe]=+0.2 and an age of 8 Gyrs. We detect no significant differences between the two populations, 
calculating median [Fc/H] and [Fe/H] t ^ ic j, = —0.8 ± 0.2. Centre panel: increasing the age of isochrones used to calculate metallicity 
for the thick disc from 8 Gyrs to fO Gyrs. An offset of ~ 0.05 dex in the average [Fe/H] of the two components is observed, with median 
[Fe/H]thick = —0.85. The dispersion remains the same as before. Right panel: Increasing [o/Fe] for the thick disc from +0.2 to +0.4. 
An offset of ~ O.f dex between the median metallicities of the populations is now observed, with [Fe/H] t ^ ic fc = —0.93. 



neighbourhood of 9:1 (|just fc Jahreifj I2010T ) , giving lxlO 6 
stars, and we use a thick disc scale length of 8.0 kpc, as 
determined above. For the velocities, we assume the thick 
disc lags behind the thin by 50 kms -1 and has the same 
velocity dispersion as the MW thick disc, a = 40 kms -1 
ijlvezic et al.ll2008h . We also generate vertical heights within 
the discs for both thin and thick populations, assuming MW 
scale heigh ts for the discs (30 pc and 1000 pc for thin and 
thick disc. Ilvezic et al1l200ct) . /-band magnitudes between 
25.0 > I > 20.3 (0.1 mags brighter than the tip of the RGB 
in M31) and angular positions within the disc. We then con- 
vert our disc frame velocities into heliocent r ic velo cities us- 
ing the HI rotation curve of IChemin et alj l|2009l '). Finally, 
we interpret this model in the same way as our data, by ro- 
tating it into the coordinate system of M31 as observed from 
the MW (with inclination and PA as discussed in § 3), and 
subtracting off the assumed disc velocity at that position by 
interpolating each of our model stars into our average disc 
velocity map (Fig. H}. 

From this model data set, we select stars as we would 
select targets to observe when designing DEIMOS masks, 
requiring them to have /-band magnitudes between 22.0 > 
/ > 20.5. We then randomly select the same number of stars 
as are observed at each field location, and make velocity his- 
tograms in the disc-lag frame for each field. We then analyse 
these distributions with the same GMM technique described 
in § 3.1, using a LRT to determine whether the distribution 
of each model field is best fit by a single thin disc compo- 
nent, or a double thin+thick component. This procedure is 
repeated 100 times, allowing us to compute the average ve- 
locities and dispersions for each component, plus sampling 
errors which we tabulate in Tabled In our final 100 samples, 
the thick disc is detected in 15 of the 21 fields on average. 
The fact that we do not see the thick disc component in all 
our model fields implies that the non-detections in our data 
are an effect of our sampling of the DEIMOS fields rather 
than the component being absent in these fields. We show 



the histograms and best fit Gaussians for three of these re- 
alizations compared to our data in Fig. 

We now assess how both the lag between components 
and the dispersions of each component evolve with radius for 
our model data set, and how accurately we can recover these 
values from our model. In Fig. 1161 we plot these values for 
our model (black circles) alongside the values we obtained 
from our data (red squares), and fit the evolution of the 
model results with linear functions as before. In the top 
panel, we show the measurement of Aw for our model fields 
as a function of projected radius. The model results show 
no evolution of Av with radius, recovering an average lag 
across all fields of 48.9+6.7 kms -1 , which is similar to the 
constant lag of 50 kms -1 implemented in our model. 

Next, we compare the evolution of the disc dispersions 
for our model with the data, shown in the central and lower 
panels of Fig. [TS] The model thin disc is best fit with a 
constant relation, giving an average lag of athin = 21.5 ± 
1.7kms _1 , very close to the input of 25.0 kms -1 . For the 
thick disc dispersion, a t hick the measured dispersion in the 
model scatters about a mean dispersion of athick = 38.8 ± 
2.4 kms -1 , with no evidence of evolution with radius. 

Finally, we can use our model to get a handle on how 
accurate our estimates of the scale length of the M31 discs 
might be. We use the same Monte Carlo (MC) technique 
above to calculate the density of stars in each component 
in our model fields 100 times, then we compute the aver- 
age density from these results. These results are shown in 
Table [6] and the errors represent the dispersion of the den- 
sities computed in each field. We then plot the densities as 
a function of radius for the thin and thick discs (shown in 
Fig. 117) 1. and fit the result with an exponential profile to 
determine the scale length. For the thin disc, we calculate 
h,- = 6.2 ± 0.8 kpc, which is consistent with our input of 6.6 
kpc, and for the thick disc we compute h r = 7.8 ± 0.9 kpc, 
consistent with our input of 8.0 kpc. The shaded regions 
indicate the la uncertainties from the fit. These results sug- 
gest that our observationally derived scale lengths for the 
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Table 5. Average kinematic properties of model fields from 100 MC realisations 



Field 
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223Mod 
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Table 6. Average densities for thin and thick discs in model fields 
from 100 MC realisations 



Field p,thin (*/arcmin) p*thick (*/arcmin) 


228Mod 


73.0±35.2 


31. 0± 19.2 
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166Mod 


14.2 ± 7.7 
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106Mod 


11.9 ± 5.7 
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224Mod 
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232Mod 


10.9 ± 6.1 


7.6 ± 4.5 


104Mod 


10.7 ± 4.9 


7.9 ± 3.8 


220Mod 


4.7 ± 1.8 


4.2 ± 1.6 


213Mod 


4.6 ± 2.6 


4.4 ± 2.4 


102Mod 


3.7 ± 1.5 


3.8 ± 1.5 


231Mod 


3.8 ± 1.9 


2.8 ± 1.3 


223Mod 


2.7 ± 1.0 


2.7 ± 0.9 


lOlMod 


2.9 ± 1.1 


2.8 ± 1.0 


222Mod 


2.1 ± 0.7 


2.1 ± 0.6 


221Mod 


1.7 ± 0.5 


1.9 ± 0.5 


50Mod 


1.3 ± 0.3 


1.2 ± 0.2 


107Mod 


1.4 ± 0.3 


1.5 ± 0.3 


wllMod 


1.4 ± 0.4 


1.5 ± 0.3 


167Mod 


0.7 ± 0.4 


0.5 ± 0.4 


148Mod 


0.6 ± 0.4 


0.5 ± 0.3 


and thick discs 


are a good 


indicator of their true i 



lengths. 

5.2 Comparison to the MW and 'edge-on' thick 
discs 

Now that we have characterised the radial profile, kinematics 
and metallicity of the thick disc in M31, we are able to com- 
pare it to the properties of other thick discs that have been 
observed in the universe. We shall begin with the most well 



studied thick disc currently known - that of our own Galaxy. 
Given that these two galaxies are relatively close to one an- 
other (separated by 785 kpc), and have similar morphologies 
(both large spiral galaxies), comparisons between the MW 
and M31 are often made. But for all their apparent similari- 
ties, thes e two galaxies a r e quit e different from one another. 
Work by lHammer et all (|2007T l has shown that the MW is 
quite different in terms of its structure and evolutionary 
history from the majority of local spiral galaxies, whereas 
M31 is actually quite "typical", so these differences are per- 
haps unsurprising. In this work, we have demonstrated that 
the scale lengths of the M31 discs are larger than those of 
the MW by a factor of ~ 2, as shown in Fig. 1101 Given 
that we derive the scale heights of the M31 discs from these 
scale lengths, this results in scale heights in M31 that are 
of order ~ 3 times as thick as those of the MW. However, 
we note that as we calculate scale heights for the M31 disc 
based on a relation determined from disc galaxies that are 
quite different in terms of their mass to both the MW and 
M31, our values may be an overestimate. The M31 discs are 
also seemingly hotter than the MW disc s, with athin M31 = 
32.0 kms -1 cf. a thi n,Mw = 20.0 kms -1 jlvezic et aljl2008h . 
and trthick.M3i = 45.7 kms" 1 cf. a th ick,MW = 40.0kms _1 
l|lvezic et al.ll2008h . This could tell us something about the 
merger history of M31. If the thick discs in both galaxies are 
formed as a result of heating by mergers, the hotter discs of 
M31 could imply that this galaxy has undergone a more 
active merger history than the MW. 

The MW thick disc is more metal poor, enriched in a 
metals and older than the thin disc. While we are unable 
to measure the age and a abundances of the M31 discs, 
we have shown that there exists an offset in the average 
metallicities of the two components of ~ 0.2 dex when mea- 
sured spectroscopically. While we do not see this offset pho- 
tometrically, this could be due to our analysis technique as 
we use isochrones of the same a abundance ([a/Fe]=+0.2) 
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Figure 15. A comparison of the data (top left panel) with 3 realisations of parsing our thin + thick disc model through the same analysis 
as our data, selecting stars from the same regions as the data. It can be seen that the model data resembles the actual data very closely, 
and that non-detections are likely an effect of sampling. 



and age (8 Gyrs) for both components. If we modify the 
a-abundance and age of these isochrones to [a/Fe]=+0.4 
and 10 Gyrs for our thick disc sample, we see an offset 
of ~ 0.2 dex. We also note that both the thin and thick 
discs in M31 appear to be more metal-poor than the MW 
discs, which have average metallicities of [Fe/H]~ —0.3 
and [Fe/Hl~ -0.6 (|Gilmore et all 120021 : lAbadi et alj|2003l ; 
ICarollo et al.l[201oh respectively, altho ugh there is a signif- 
icant metallicity spread in both discs. ICarollo et al.l (|20ld ) 
also demonstrated evidence of a secondary, more metal poor 
thick disc in the MW, whose metallicities span the range 
-1.8 < [Fe/H] < -0.8, peaking at [Fe/H]=-1.3. This com- 
ponent also appears to be hotter than the traditional MW 
thick disc component, with cr z =44±3kms~ 1 , very similar to 
what we observe in M31. 

In §4.1.3, we inferred scale heights for the thin and thick 
discs of M31 of h z = 1.1 ± 0.2 kpc and h z = 2.8 ± 0.6 kpc 
respectively, using a sample of 34 galaxies with thick discs 
measured by YD06 to determine a relationship between scale 



length and scale height of a stellar disc. As we noted in 
§4.1.3, a comparison with the YD06 sample might not be 
desirable, as these galaxies are typically much less massive 
than M31, and selected to be bulgeless. A more appropriate 
comparison would be the MW analogue, NGC 891, an edge- 
on gala xy that was r ecent ly the subject of a structural anal- 
ysis by llbata et all (|2009l ) using HST/ACS imaging. They 
detected the presence of a thick disc component in the galaxy 
and were able to measure both a scale length and height for 
this component of h r = 4.8±0.1 kpc and zo — 1.44±0.03 kpc, 
compared with h r = 4.2 ± 0.01 kpc and zo = 0.57 ± 0.01 kpc 
for the thin disc component in this galaxy. This gives a ratio 
of ~ 1.1 between the scale lengths and ~ 2.5 for the scale 
heights of these components, which is identical to what we 
observe in M31. To illustrate this, we overplot these values 
for NGC 891 in Fig. [9] as a green circle. 

In lYoachim fc Dalcantonl l|2008bh the authors present 
kinematics of the thin and thick discs of 9 of their initial 
sample of 34 galaxies, obtained using the GMOS spectro- 
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Figure 16. This figure compares the results from our data with 
results from our model analysed in the same way. In all cases, 
data is represented by filled red squares and dot-dashed lines, 
and the model results are shown as filled black circles and solid 
lines. Top panel: The difference in velocity, Av, between the 
thin disc and thick component of data and model as a function of 
projected radius. The model lag is consistent with no evolution 
with radius, and shows an average lag of 48.9 kms , very close 
to our input lag of 50 kms -1 . Middle panel: Dispersion, dtuini 
of the thin disc is plotted for both data and model as a function of 
projected radius. The model thin disc is best fit with an average 
dispersion of 21.5 kms -1 , very close to the input of 25.0 kms - . 
Lower panel Results for both data and model for the dispersion 
of the thick disc (pthiak) as a function of radius. For our model, 
the thick disc is consistent with no evolution with radius, unlike 
our data, and has an average dispersion of a v = 38.8 kms -1 , 
which recovers our input dispersion of 40 kms - relatively well. 



graph on Gemini. To measure velocities and dispersions in 
both thin and thick components, they placed slits in posi- 
tions corresponding to the midplane of the galaxy to mea- 
sure the thin disc properties, and above the midplane where 
the contribution from the thin disc was thought to be neg- 
ligible. As their typical velocity resolution was 60 kms -1 , 
they were unable to draw robust conclusions on the veloc- 
ity dispersions of these components, but they were able to 
measure velocity rotation curves for each component, and 
found a wide variety of behaviour amongst their thick disc 
components, with discs which lagged behind the thin disc 
by only ~ 5 kms - , discs that show no evidence of rota- 
tion and one case where the thick disc is counter-rotating 
with respect to the thin disc. The average lag between the 
thin and thick components of Av = 46.0 kms -1 we see in 
the M31 system is larger than the majority that they ob- 
serve. We note that the galaxies in their sample were typi- 
cally of much lower mass than M31 (V circ < 150 kms -1 cf 
Vdrc ~ 230 kms -1 ). In the most massive of their sample 
(which are still less massive than M31), they do not detect 
a lag in the thick disc kinematics at all, and they attribute 
this to contrast issues. Their sample were also selected to be 
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Figure 17. Results from MC recovery of the scale lengths in 
our thin+thick disc model for our input thin (top panel) and 
thick (bottom panel) discs. The error bars on individual points 
represent the dispersion of calculated densities in the MC analysis, 
while the shaded regions represent the la uncertainties from the 
weighted least-squares fit. We recover a scale length for the thin 
disc of h r =6.2±0.8 kpc and h r = 7.8±0.9 for the thick disc, which 
are consistent with our input values. 



"bulgeless", unlike M31 which has a significant bulge, and 
so a direct comparison may not be advisable. Owing to the 
wide range of kinematic behaviour exhibited in their sam- 
ple, they conclude that the dominant formation process of 
th ick discs is via minor mergers a nd accretions of satellites. 
In lYoachim fc Dalcantonl (|2008au . they use Lick indices to 
measure ages and metallicities in 9 low mass galaxies with 
thick disc components. While we measure an offset of 0.2 dex 
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in the metallicities of the M31 thick and thin discs, they were 
unable to measure any such offset in their sample, though 
this could be a result of the insensitivity of Lick indices to 
such differences at low metallicity. They do find that the 
thick discs are host to older stellar populations than the 
thin disc, however with our current data set, we are unable 
to comment on the ages of stars in the M31 discs. 

5.3 Possible formation scenarios 

In this section, we discuss the various formation scenarios 
mentioned in § 1. Owing to our inability to measure ages 
and vertical dispersions in M31, we are not able to confirm 
or reject any of these formation mechanisms at present, so 
we discuss additional constraints for these models that could 
help to rule out or confirm each scenario with further data 
and analysis. 



5.3.1 Heating by minor mergers 

Numerous studies have identified that impacts and merg- 
ers of satellites with masses less than a third of their hosts 
can kinematically heat the thin stellar disc, puffing it out 



Robin et all Il996l; Iwalker et all Il996; 
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; Chen et al. 2001; Sales et al.I:2009; 
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2009). M31 is known to have recently undergone at least one 
significant minor merger event, resulting in t he GSS tidal 
stream. In recent work bv lPurcell et alj |2010T l. the authors 
model the heating of the stellar discs by minor mergers and 
trace disc stars ejected into the stellar halo by these simu- 
lated events. In addition to the stars ejected into the halo, 
they observe a concomitant increase in the number of stars 
located in the kinematic regime of the thick disc, contribut- 
ing ~ 10 — 20% of the total stellar density along the major 
axis, similar to what we observe in M31. They also find that 
their simulated planar infall produces two-component sys- 
tems with scale heights (z t Mn ~ 1 kpc and z thick ~ 3 kpc), 
consistent with our measurement s for M31. The simi larities 
between our findings and those of lPurcell et al] (|2010h could 
suggest that thin disc stars heated by the merging event that 
created the GSS may contribute some non-trivial fraction of 
stars to the thick disc. 

According to the simulations of iKazantzidis et all 

(2009), thick discs produced in this vein imprint a number of 
dynamical signatures on both the kinematic and structural 
properties of the galaxy. These include considerable thicken- 
ing and heating at all radii, prominent flaring, particularly 
in the outskirts of the disc (beyond 3 scale lengths), sur- 
face density excesses at large radii, radial anisotropies and 
substantial tilting of the disc. As M31 is not edge on, we 
are unable to comment on the evolution of the height of the 
thick disc with radius, and so we cannot use this as a mea- 
sure of flaring in the outer regions of the disc. However, one 
might expect that if there was a substantial flaring beyond 3 
disc scale lengths (~ 24 kpc), that this may be reflected by 
an increase in the velocity dispersions of both thin and thick 
disc components. Our results for evolution in the thin and 
thick disc dispersions remain inconclusive, and so it is pos- 
sible that such flaring may exist. At present, we possess few 
fields between R~ 32 and 39.6 kpc, so populating this re- 
gion with kinematics, as well as additional fields further out, 



may further enlighten us to any potential flaring. Another 
test of this formation scenario would be to include fields 
from both the minor axis and NE portion of M31 to test for 
any radial anisotropy, assuming one can reliably disentan- 
gle contamination from foreground and substructure from 
the signatures of the discs. The work of lSales et al] (|2009l ) 
also tells us that thick discs that are produced as a result of 
heating present structures with low orbital eccentricity. 



5.3.2 Accretion of satellite on a coplanar orbit 

Numerical simu l ations by lAbadi et al] (|2003l ) and 
iPenarrubia et aL (2006) show that an old, thick disc 
of stars could form via the accretion of stars from satellite 
galaxies on an approximately coplanar orbit with its host. 
Such discs are similar in radial extent and contain older 
stellar populations when compared to the thin disc. The 
thick disc we find in M31 is consistent with this model in 
so far as the radial extents of both discs (5.9-7.3 kpc and 
8.0 kpc) are comparable with one another. They also argue 
that the mass and luminosity of the progenitor satellite 
can be inferred from the metallicity of the component. We 
deduce [Fe/H] tMcfc = -1.0 ± 0.1 for M31, which would 
correspond to a satellite of My ~ —15 (L v ~ 9 x 10 7 Lq), 
similar to the M31 d warf elliptical, NGC 147 (M v = -15.1, 
Ivan den Bergh 1999). However, given the mass we calculate 
for the thick disc in §4.1.4 of 2-4x 10 10 M , it seems very 
unlikely that the thick disc of M31 could have been formed 

from such a satellite. 

Results of the simulations of ISales et al] (I2009D show 
that stars accreted into a thick disc from satellites on copla- 
nar orbits exhibit high eccentricity orbits. Our present data 
set does not allow us to probe the eccentricity of the orbits 
within the thick disc at this time. With a larger data set, 
we could perhaps see the effects of orbital differences in the 
form of structural asymmetries. 



5.3.3 Radial migration and internal heating 

The scattering of stars by spiral structure and molecular 
clouds has long been proposed as a method of heating 
the stellar disc, moving stars out onto m o re eccentr i c and 
inclined orbits (ISellwood fc Binnevl |2002|; iHavwoodl 120081 : 
iRoskar et ailbOOSl; ISchonrich fc Binnevl l2009a|) . and it has 
been argued in ISchonrich fc Binnevl (|2009a| ]bh that these 
processes naturally produce an old, a-enhanced thick disc, 
whose properties are consistent with those observed in the 
MW. These models also show wide MDFs and an increase 
in the scatter of the Age-Meta l licity relation. This is also 
demonstrated in lOuillen et al] (|2009l ) , where they investi- 
gate radial mixing induced by an orbiting subhalo. Again 
they find evidence of wide MDFs in both the thin and thick 
discs. With deeper photometry that allowed us to reach 
the MSTOs of the two discs we could derive the average 
ages of these components, and high resolution spectroscopy 
(R~ 15, 000) of M31 thick disc stars that would allow us to 
determine accurate abundances from unblended Fe lines for 
individual stars, we could comment more robustly on the 
likelihood of such a formation scenario. 
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5.3.4 Thick disc forms thick 

iKroupal (|2002h posited that thick discs could be formed as 
a result of vigorous star formation in massive star clusters 
(~ 10 5 — 10 6 Mq) during the period of assembly of the stel- 
lar disc. If this is true, a number of these massive clusters 
may have survived to the present day, an d woul d possess 
large vertical velocity dispersions. IKroupal (2002) suggests 
that these clusters could be the metal-rich globular cluster 
system in the MW. Once again, owing to the inclination of 
M31, we are unable to measure the vertical dispersions of its 
metal-rich globular cluster system, and can therefore neither 
confirm nor reject this formation model. 



6 CONCLUSIONS 

Using the DEIMOS spectrograph on the Keck II telescope, 
we have identified a statistically significant population of 
stars in M31 that lags behind the thin and extended discs 
by 46.0±3.9 kms^ 1 . Comparing this with a model of a 
thin+thick disc system with the same distance and incli- 
nation as M31 shows this component to be consistent with 
a thick disc component. Analysing its kinematics, we find 
it to be hotter than the thin disc, with average dispersion 
o'thick = 50.8 ± 1.9kms _1 cf. a thin = 35.7 ± l.Okms , 
larger than the dispersions observed in the MW discs. From 
composite spectra for each component, constructed from 
highly probable thin and thick disc stars (selected using 
stringent Gaussian cuts) we measure a metallicity offset of 
~ 0.3 dex between the two disc, with the thick disc being 
metal-poor than the thin disc ([Fe/HL, . , = —1.0 ± 0.1 cf. 
[Fe/H] tMn = -0.7±0.05). The fact that'this metallicity off- 
set is not observed when analysing the thin and thick disc 
RGB stars with isochrones of identical age and a-abundance 
suggests that the two populations differ in these properties, 
with the thick disc likely being older and more enriched in 
a elements. 

We measure scale lengths for both thin and thick discs, 
finding h r = 8.0 ± 1.2 kpc for the thick disc, and h r = 7.3 ± 
1.1 kpc for the thin disc, comparable to previous estimates. 
Using the data of YD06 we infer scale heights for both discs 
at zo = 2.8 ± 0.6 kpc and zq = 1.1 ± 0.2 kpc for thick and 
thin discs respectively. These values are of order 2-3 times 
larger than those measured in the MW, perhaps suggesting 
that M31 has undergone more heating than our Galaxy. 

By measuring the ratio of the densities of both discs, we 
are able to estimate a mass range for the thick disc compo- 
nent of 2.4x 10 10 M Q < M» )t hick < 4.1 x 10 10 M Q . This value 
provides a useful constraint on possible formation mecha- 
nisms, as any proposed method for forming a thick disc must 
be able to heat (or deposit) at least this amount of material. 

Owing to current limitations within our data set, we are 
not able to distinguish between the different thick disc for- 
mation mechanisms. However, with further analysis of this 
component using our complete kinematic sample (including 
regions from the minor axis and NE of M31) and spectro- 
scopic follow up of fields where this component is strongly 
observed, we will be able to better understand the chemistry 
of this component and distinguish between various forma- 
tion mechanisms. 
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